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Abstract 

Possible mechanisms for gluon confinement in Yang-Mills theory in covariant gauges 
are discussed briefly. Implications for the infrared behavior of the Green's functions 
are outlined. Results for the gluon and Faddeev-Popov ghost propagators from Dyson- 
Schwinger studies and lattice calculations are presented. The importance of the ghost- 
gluon vertex in these studies is discussed. Numerical results for a non-perturbatively 
constructed ghost-gluon vertex in four and three space-time dimensions are presented. 



1 Confinement in Covariant Gauges 

The strong interactions governing hadron and ultimately nuclear dynamics are described 
by quantum chromodynamics, the theory of quarks and gluons. This theory has several 
genuine non-perturbative features, such as e.g. dynamical chiral symmetry breaking. The 
most prominent phenomenon is, however, confinement. Strong evidence exists that it is 
generated in the Yang- Mills subsector, i.e. by the gluons. Naively, a Yang-Mills theory is 
defined in terms of the generating functional P] 



z[j] = j vAi exp j d'^x c + j d'^x 



r ^ rpa rpa 

F;, = d,Al-d.Al + g,r"^AlA 
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Here we have already chosen a description in Euchdean space-time, which is used throughout 
this work. Aj^(a;) is the gluon field, F^^{x) is the field strength tensor, is the gauge 
coupling constant in d dimensions, and /'^'"^ are the structure constants of the gauge group. 
The measure of the functional integral extends over the complete configuration space of 
the gluon field. 

However, the functional integral (0) is ill-defined. The Lagrangian is invariant under 
local gauge transformations and therefore the complete configuration space contains infinite 
sets of mutually equivalent field configurations. They are said to lie on the same gauge orbit. 
To avoid this over-counting in the functional integral, it is necessary to select only one copy 
on each gauge orbit. One tries to achieve this by gauge-fixing ilj. Implementing the linear 
covariant gauge fixing condition d^A^^ = leads to the generating functional 

Z[J] = jvA^M exp j d'x iC + C,i) + J d'xj;A';}j , 
M = det{-d,Df), (2) 



= aj'^' + g^r'^A 



where M is the Faddeev-Popov determinant, and £gf the gauge-fixing term. The determinant 
is then rewritten as a functional integral over scalar Grassmann fields, the Faddeev-Popov 
ghosts. The generating functional then reads 



(3) 



The ghost and anti-ghost fields, c°-{x) and c"-{x), respectively, encode part of the quantum 
fiuctuations of the gluon field. Their appearance is a generic feature of covariant gauges. As 
will be detailed below, this procedure is in general not sufficient to completely specify the 
gauge beyond the use in perturbation theory. 

One of the most intriguing features of QCD is confinement. Mechanisms for gluon con- 
finement have been related to the infrared behaviour of propagators, see e.g. ref. We 
will focus here on aspects which have been substantiated recently, cf. section |21 

First we discuss a sufficient but not necessary condition for confinement. If a propagator 
D vanishes sufficiently fast in the infrared, 

lim D{p^) = , (4) 

then the "particle" represented by D cannot have a positive semi-definite spectral norm. 
Therefore it violates the Osterwalder-Schrader axiom of refiection positivity and the 
propagator D cannot have a Kallen-Lehmann representation. Thus this "particle" is not 
part of the physical spectrum, and is hence confined 

For a massless excitation such as the gluon this confinement criterion is rather intu- 
itive. Then, the condition (j^ is just the statement that the on-shell propagator vanishes. 
Therefore, such a "particle" does not propagate, i.e. it is confined. 

The criterion (0} is related to a scenario which was introduced by Kugo and Ojima 
In perturbation theory, ghosts and longitudinal gluons are not part of the physical spectrum, 
since they form doublets with respect to BRST charge. As these doublets are conjugate pairs 
by ghost number, this is called the quartet mechanism. The conjecture of Kugo and Ojima 
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states that also transverse gluons form together with ghost-gluon bound states a quartet, 
and are thus not part of the physical spectrum. This does not imply that transverse gluons 
do not have asymptotic states. However, these states are not in the physical subspace and 
therefore do not contribute to ^-matrix elements. 

This scenario is based on a proof which requires three preconditions. The first is an 
unbroken BRST charge. Whether this is the case is not known. The second is the failure of 
cluster decomposition in the complete state space. This requires the existence of a massless 
excitation in the unphysical subsector. However, at the same time a mass gap in the physical 
subsector is required. The third precondition is an unbroken global color charge. This 
condition can be cast in Landau gauge into a condition for the ghost propagator Dq [HI 



As will be seen below this scenario also entails condition for the gluon propagator. 

Another scenario has been proposed by Zwanziger [7j based on work by Gribov [S]. Only 
one gauge copy on each gauge orbit should be included in the functional integral constituting 
the generating functional. However, no local gauge condition is known which is able to specify 
the gauge completely also for large gauge field fluctuations. This is known as the Gribov 
problem [S]. The origin of this problem actually provides a possible handle to understand 
confinement. Therefore a closer look at the properties of the field configuration space is 
necessary. 

It is possible to partition configuration space by the zeros of the Faddeev-Popov deter- 
minant (j2I). Concentric regions are found and each gauge orbit intersects each region at 
least once jH|. The region containing the origin, and thus relevant for perturbation theory, is 
known as the first Gribov horizon. It can be shown to be compact and convex. It has been 
found, that this region also contains gauge copies, see e.g. PJ, and it is necessary to further 
restrict to the so-called fundamental modular region. This region is specified by minimizing 
an appropriate functional along each gauge orbit ^liij. The fundamental modular region is 
also convex and compact around the origin [7j. 

The latter properties lead to an interesting observation. As configuration space is infinite- 
dimensional, any volume is dominated by its boundary. It has therefore been argued that 
only the common boundary of the first Gribov horizon and the fundamental modular region 
contribute to any correlation function built from a finite number of field operators, especially 
Green's functions [7j. Therefore, restricting to the first Gribov horizon, would be sufficient 
to calculate Green's functions. 

This has a further important consequence [7j. The Faddeev-Popov determinant Q van- 
ishes by definition on this boundary. Therefore, the condition (0), termed in this context 
the Zwanziger horizon condition, is fulfilled as the determinant is directly linked to the in- 
verse self-energy of the ghost. This also entails that the action is dominated by the gauge 
fixing term alone. Thus, in the extreme infrared all Green's functions are determined by 
contributions involving ghosts. This leads to the notion of ghost dominance. 

Assuming the horizon condition (0), it can be shown that an infrared tree- level ghost- 
gluon vertex complies with the renormalization group in leading order. This leads to the 
hypothesis of a bare ghost-gluon vertex also in the infrared. It has then been shown that, 
for the Landau gauge gluon propagator, the condition (jH) follows and gluons are confined 




(5) 
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Figure 1: The Dyson- Schwinger equations for the 2-point functions of Yang-Mills theory. 
Dotted lines are ghosts and wiggly lines are gluons. Lines with a dot are full propagators, 
vertices with a small black dot are bare vertices and white circles are full vertices. The 
equation for the gluon has already been truncated to the one- loop level, see text. 

Whether the identical infrared behavior of the ghost and gluon propagators within the 
Kugo-Ojima and the Zwanziger-Gribov scenarios, respectively, is more than a formal coinci- 
dence is not known. However, a ghost propagator with a stronger infrared divergence than 
a massless pole is also intuitively seen to be able to mediate confinement. After Fourier 
transformation such a singularity corresponds to a correlation being more long-range than 
that of a photon, and thus mediates a stronger force than a Coulomb force. 



2 Propagators in Yang-Mills Theory 

In both confinement scenarios described above the infrared behavior of propagators plays a 
central role. As infrared singularities are anticipated, cf. eq. (0), a non-perturbative method 
formulated in space-time continuum is desirable. Such a formulation is provided by the 
Dyson-Schwinger equations (DSEs), determining the Green's function of a theory, see e.g. 
PlEj. A graphical representation of the equations for the propagators is shown in fig. E 

The Landau gauge will be used throughout, as it is advantageous in the study of con- 
finement for various reasons. First, the Landau gauge is a fixed point of the renormalization 
group. Second, the conditions for the Kugo-Ojima and Zwanziger-Gribov scenario take the 
simple forms (jU and (0) in Landau gauge. Last but not least, it is less singular than 
other gauges. This manifests itself in the bareness of the ghost-gluon vertex for a vanishing 
incoming ghost momentum ^Tj 

]imTf'^{k;q,p)=zg,r''^q,. (6) 
p— >o ^ 

Here p, q, and k are the momenta of the incoming ghost, the outgoing ghost, and the gluon, 
respectively. 

Employing the ghost DSE and quite general assumptions, it can then be shown [12j that 
the ghost propagator necessarily fulfills condition (0). It is not possible in this generality to 
show that the gluon propagator satisfies (jlj). This requires a solution of the DSEs. However, 
it is not possible to solve these equations as they form an infinite set of coupled non-linear 
integral equations. Therefore, it is necessary to truncate the system. 

During the last years, the truncation scheme depicted in fig. ^ has been established JH] ■ 
This scheme has been successfully applied to Yang- Mills theory in four dimensions [Tij . 
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three dimensions [T3] and at finite temperature ^HI as well as to QCD including dynamical 
quarks [T7]. In this scheme, all equations beyond those of the propagators are neglected, as 
well as genuine two-loop contributions in the equations for the propagators. Furthermore, 
a perturbative color structure of the Green's functions is assumed. It then only remains 
to specify the two remaining full vertices, the ghost-gluon and the 3-gluon vertex. The 3- 
gluon vertex is chosen such that leading-order resummed perturbation theory is reproduced 
in the ultraviolet. Due to the property (jB]) the ghost-gluon vertex is usually taken to be 
bare. The results for the ghost and gluon propagators are shown in figs. El and El compared 
to lattice calculations. Both methods agree remarkably well. The largest deviations are 
found at intermediate momenta where the truncation scheme is expected to be least reliable. 
Although this scheme imposes several problems, especially concerning gauge invariance, in a 
systematic analysis only small quantitative effects due to truncation errors have been found, 
for a detailed discussion of this issue see e.g. refs. jTHHnj. 




Figure 2: The ghost dressing function from DSEs in four dimensions compared to lattice 
results ^Hl (left) and from DSEs in three dimensions P3] (right), is the dimensionful 
coupling constant in three-dimensional Yang-Mills theory. 

The gluon propagator obeys condition (0]), and the gluon is thus confined. As the 
ghost propagator also satisfies the Zwanziger horizon condition (jSJ, the Kugo-Ojima and 
the Zwanziger-Gribov scenarios are both fulfilled. These findings have also been confirmed 
recently by exact renormalization group methods [22] • Furthermore, such propagators are 
in sharp contrast to the older concept of infrared slavery, which expected a gluon dressing 
function diverging as l/p^ [22 • These earlier attempts neglected the ghost contribution 
pij . This so-called Mandelstam approximation leads to results contradicting current lattice 
results and can now be dismissed. 

In the infrared, ghost dominance is found. Thus, besides the full gluon propagator, only 
diagrams with at least one ghost line in fig. ^ contribute. This is in accordance with the 
Zwanziger-Gribov scenario. The single crucial point in this derivation is then the structure 
of the ghost-gluon vertex j2Sj- It is therefore of great importance to understand it. The 
infrared behavior of the propagators has been studied using several ansatze, motivated by 
the Slavnov- Taylor identity (STIs) for the ghost-gluon vertex. The results are found to 
depend only quantitatively on the input vertex [25j. However, a more extensive analysis is 
highly desirable and will be presented below. 
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Figure 3: The gluon propagator from DSEs in four dimensions (HI compared to lattice data 
[TU] (left) and from DSEs in three dimensions jl5j compared to corresponding lattice data 
[201121] (right). Note that the lattice data at low momenta are likely to be strongly influenced 
by finite volume effects. 



3 Ghost-Gluon Vertex 

The most general tensor structure of the ghost-gluon vertex in Landau gauge is given by 

Tf{k; q,p) = ig, {q^ {^'^ + A'^'^k; q,p)) + k^B'^'^k- q,p)) , (7) 

where A°''"^ and B"''"^ are scalar functions. The tree-level contribution is left explicit. As 
there is no indication for a color structure different from perturbation theory [2. , the struc- 
ture constant is factored out, A"'^'^ =: f'^^'^A and B"''"^ =: f°''"^B. Due to the transversality 
of the gluon propagator, the on-shell ghost-gluon vertex is transverse with respect to the 
gluon momentum. Therefore, the function B is only relevant when not contracting with 
a transverse projector, e.g. in the corresponding STL Hence, only the function A will be 
discussed here. A detailed account on B is given elsewhere P"6J. 

In principle, the ghost-gluon vertex should be determined self-consistently by solving the 
DSEs for the propagators and the vertex simultaneously. A truncated DSE for the vertex 
neglecting four-point functions is shown in fig. EJ The completely self-consistent solution 
of such an equation is of significant technical complexity and has not been achieved yet. 
However, due to the hypothesis of Zwanziger, it is most likely sufficient to perform a semi- 




Figure 4: The DSE for the ghost-gluon vertex. Dotted lines are ghosts and wiggly lines 
are gluons. Lines with a dot are full propagators, vertices with a small black dot are bare 
vertices and white circles are full vertices. Contributions from the ghost-gluon scattering 
kernel have been removed, see text. 
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Ghost-gluon vertex in 4d 
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Figure 5: The normalized transverse part (jH)) of the ghost-gluon vertex in four dimensions 
in various kinematical regions. 



perturbative calculation, i.e. to do one iteration step in the ghost-gluon vertex DSE. If the 
hypothesis is correct then the resulting vertex should not significantly deviate from the input 
tree-level vertex. 

Hence, the evaluation scheme neglects the ghost-gluon scattering kernel in the ghost- 
gluon vertex DSE. Further, as this is the first iteration step, all vertices are left bare and 
only the propagators are dressed with the results presented in the previous section. The 
results are displayed in fig. El for four dimensions and in fig. IHlfor three dimensions. Shown 
is the transverse part of the ghost-gluon vertex, 1 + A, defined by 

- ^) <i^p) = ^'^(1 + Mk; q,p)) (8) 

where A = q'^k^ — {q ■ k)"^ is a Gram determinant. 

Several observations can be made First, although not shown in the plots, the 

results comply with the infrared limit (jHI) and thus satisfy the respective STI in this specific 
kinematical limit. Second, they have the proper ghost-antighost symmetry = 
A{k;p,q). Therefore, they fulfill all required conditions. In addition, in three dimensions, 
asymptotic freedom is manifest and the vertex reduces to the bare one in the ultraviolet. 
However, in four dimensions the anomalous dimensions is incorrect, as the vertex is not 
resummed yet. 

As can be seen, the result is a vertex which is very close to the tree-level one at all 
momenta, especially in the infrared. This does not change in other kinematical conditions 
and thus confirms Zwanziger's hypothesis in an impressive manner. The origin for 
the small deviations from the exact bare vertex in the infrared may be twofold. On the 
one hand, certainly iteration to self-consistency and truncation artifacts change the results 
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Ghost-gluon vertex in 3d 
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Figure 6: The normalized transverse part (jH)) of the ghost-gluon vertex in three dimensions 
in various kinematical regions. 



quantitatively. As the ghost-gluon vertex DSE is already quite well fulfilled, this will most 
likely not lead to a qualitative change. On the other hand, in Zwanziger's hypothesis only 
agreement with leading-order renormalization group has been required. Subleading effects 
may generate a deviation of the order observed. 

In addition, the vertex is nearly tree-level for all momenta. This supports the assumption 
of a bare ghost-gluon vertex in the truncation scheme for the propagators. It especially 
provides additional evidence for ghost dominance and therefore for the Zwanziger-Gribov 
confinement mechanism. 

Similar results have also been reported recently in lattice calculations in four dimensions 
P7j . They confirm the results found independently here and indicate also a nearly tree-level 
ghost-gluon vertex. 

To assess the error of the approximations adopted, the full ghost-gluon vertices in the 
ansatz shown in fig. E] have been replaced by the vertices motivated by the corresponding 
STI It has been found [20] that also these vertices are mapped to a nearly tree-level 
like vertex by the truncated DSE, as long as no singularity in an external momentum is 
introduced. (As there is no iteration in the current ansatz, these divergencies would survive 
the calculation.) Moreover, the results have been found to be slightly more sensitive to 
the propagators than to the input vertices. All these variations are quantitatively small. 
This indicates that the solution of the ghost-gluon vertex DSE is a vertex which is nearly 
tree-level. 

It is probable that the infiuence of quarks, cf. ref. [T7j, or of an adjoint Higgs field, cf. 
ref. ^21; on these findings is negligible. Neither couple directly to the ghost and therefore 
their presence does not give rise to new terms in the ghost-gluon vertex DSE. Hence, they 
can only contribute indirectly via their effect on the Yang-Mills propagators, which are very 
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robust against coupling a small number of quark flavours J7], or scalars [T3] . 

The ghost-gluon vertex may also be interesting with respect to extracting the running 
coupling constant [2ZjJ. This issue has to be investigated further. 



4 Conclusions 

In Landau-gauge Yang-Mills theory we have calculated the ghost-gluon vertex employing 
non-perturbative propagators. The results show quite conclusively that deviations from 
the tree-level value are very small, even in the infrared. On the one hand, this supports 
the truncation scheme used to calculate the propagators of the Yang-Mills theory. The 
results found in these calculations give strong evidence for a confinement mechanism of the 
Kugo-Ojima or Zwanziger-Gribov type. Especially, no sign of an infrared enhanced gluon 
propagator in Landau gauge is found. More importantly, the results on the ghost-gluon 
vertex alone, provide independent support for the Zwanziger-Gribov scenario as they fulfill 
Zwanziger's hypothesis of a bare ghost-gluon vertex in the infrared. The results presented 
here are thus a further indication for a confinement mechanism driven by the dynamics on 
or near the Gribov horizon. This coincides with recent findings in Coulomb gauge j2Hl- AH 
these evidences combined promise a new era of quantitative confinement physics. 
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